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Abstract. We investigate a class of random graph ensembles based on the Feynman graphs of 
multidimensional integrals, representing statistical-mechanical partition functions. We show that 
the resulting ensembles of random graphs strongly resemble those defined in random graphs with 
hidden color (CDRG), generalizing the known relation of the Feynman graphs of simple one- 
dimensional integrals to random graphs with a given degree distribution. 

INTRODUCTION 

For a long time, the concept of random graphs was synomymous with the classic ER- 
model, which has been extensively studied ever since it was introduced in the sixties 
[1, 2, 3]. In the last decades, however, driven by the steadily increasing availability 
of data on large real-world networks (e.g. various social, biological, and information- 
technological networks), a large variety of more general models have been considered; 
see e.g. [4] for a fairly up-to-date review. 

One of the more commonly studied classes of models is random graphs with a given 
degree distribution, or the configuration model [5, 6, 7], here to be referred to as DRG 
(for "degree-driven random graphs"). There, one considers a maximally random graph 
with /V vertices, restricted such that its degree sequence, {m,-, i = 1 . . .N}, is consistent 
with a specified degree distribution, {p m }. Once the degrees of individual vertices are 
fixed, edges are formed by means of a random pairing of the entire set of stubs (half- 
edges), assuming an even total stub-count. This leads in general to a multigraph, with 
multiple edges and/or self-connections. The resulting graph ensembles are known to lack 
non-trivial edge correlations. 

Recently, extensions of these models based on adding so called hidden color, have 
been studied. In a colored extension of the ER-model, Inhomogeneous Random Graphs 
(IRG), vertices were assumed to carry color variables, allowed to affect edge formation, 
which was shown to result in a class of models allowing for non-trivial edge correlations 
[8], while degree distributions were limited to mixtures of Poissonians. 

Similarly, colored extensions of the DRG class of models have been considered, where 
not vertices but half-edges (stubs) were considered to carry a color variable, represented 
e.g. by an integer variable c E [1...K], allowed to affect edge probabilities. The degree 
of a vertex amounts to its stub count; since stubs are colored it is natural to consider the 
concept of a colored degree, in the form of a vector m = {m c , c = 1 . ..K}, where K is 
the number of possible colors, and m c is the number of stubs of a definite color c. The 
conventional degree m then is obtained simply by summing up the components of its 
colored degree, m = Ylc=\ m c- Thus, in the CDRG (for Colored DRG) class of models 



[9, 10, 11], random graphs with a given colored degree distribution are considered, 
where a sequence of colored degrees is chosen consistently with a specified distribution. 
The color is allowed to affect the edge distribution by introducing a color-dependent 
bias in the random pairing of stubs. The color can either be considered as an observable 
characteristic, or as unvisible once it has done its job (hidden color). 

In this article, we will study models of random graphs based on the Feynman graphs 
associated with a class of simple statistical mechanics models, where they represent the 
deviations of certain integrals from a Gaussian approximation [12]. Relations between 
DRG models and the Feynman graphs of certain simple univariate integrals are known, 
and have been explored e.g. in refs. [13, 14]. 

We will show that such relations are not unique for univariate integrals; indeed, 
any integral of whatever dimension, representing the partition function of a statistical- 
mechanical system from a large class of models, will define an ensemble of random 
graphs related to an associated CDRG model, and this is the main result of this article. 

One thing to be gained from such relations is an alternative, and perhaps more natural, 
way to define the associated DRG or CDRG random graph ensembles, in terms of 
random Feynman Graphs. 

In addition, as always when relations are found between different areas of science, one 
may expect conceptual gains, where phenomena in one area can be simply understood 
in terms of well-known phenomena in the other. Thus, the relations raise questions, e.g., 
as to how critical parameter values in the random graph models are mirrored in terms of 
phase transitions in the associated statistical physics models. 

RANDOM FEYNMAN GRAPHS - UNIVARIATE CASE 

In field theory and statistical physics, Feynman graphs provide a standard method for 
organizing the perturbative expansion of an integral Z representing some field theoret- 
ical or statistical physics model, around an approximate value as given by a Gaussian 
approximation to the integrand (see e.g. [12]). 

Sometimes this approximation is based on a saddle point, an extremal point of the 
integrand, with the Gaussian chosen to approximate the integrand well in the neighbor- 
hood of the saddle point. However, any decomposition of the integrand into the product 
of a Gaussian factor and a non-trivial remainder will do. 

In the spirit of refs. [13, 14], we will consider the possibility of defining ensembles 
of random graphs - random Feynman graphs - based on the Feynman graphs for such 
a system, with the statistical weight of each distinct graph taken to be proportional to 
the corresponding term in the perturbative expansion of Z. Much of the contents of this 
section concerns known results [13, 14], included to make the article reasonably self- 
contained. 



Feynman Graphs for a Simple Integral 



In the simplest case, the integral would be over a single real variable, x, say, with an 
integrand in the form of the exponential of an action S(x), that can be written as the sum 
of a simple quadratic term Sq(x) that defines the unperturbed action, and a non-trivial 
perturbation V(x), typically a more complicated function. Thus, consider the (formal) 
integral 

Z = [e s ^dx = Je s ^ +V ^dx =[e- h * 2 / 2 e v ^dx (1) 

where we restrict ourselves to b > 0, and to a perturbation V with only non-negative 
Taylor coefficients, 



x n 



V(x) = V v M — , with v n > 0, for n = 0, 1 . . . (2) 

n nl 

The Feynman graphs for such an integral represent different terms in the perturbative 
expansion of Z around the Gaussian approximation, Zq = J e s °^dx = J e~ bx / 2 dx 

To see this, first note that considering V(x) as a perturbation means that it is to 
be considered as small in some sense, and so it is natural to write Z as the sum 
of a sequence of contributions, obtained by expanding the non-trivial factor e v(x ) = 
Z n V(x) n /n\, yielding 



withZ^= [e- bx2 / 2 V(xfdx. (3) 

TV ' J 



Now, each of the N factors of V(x) is to be associated with a distinct vertex, labelled 
i= \ ,...,N. Next, by Taylor-expanding each factor of V(x), Z/v decomposes into a sum 
of contributions, 

Z N = Z {m,} 5 with (4) 



where the sequence {m,, i = 0, . . . ,,/V} represents the chosen terms {v m( jc m '/m ; !} in the 
respective Taylor expansions. 

Now each factor x m in the integrand represents a vertex of degree m,-, with each single 
factor of x representing a stub. Each term is now easy to evaluate. Denoting by M the 
total degree M = Y*i m u we have / e - bx2 / 2 x M dx = e M b- M / 2 {M-\)\\, where Em is a 
0/1-variable enforcing the restriction to even M; 

£ ^i±^. (6, 

The resulting factor of (Af — 1)!! has an obvious combinatorial interpretion as the total 
number of distinct pairings of the M factors of x. In each pairing, each paired couple of 
x-factors can be associated with an edge between the corresponding vertices i, j. 



By considering the m, distinct factors of x associated with the stubs of a vertex as 
indistinguishable, pairings can be naturally grouped into equivalence classes. Each such 
class of equivalent pairings can be associated with a distinct graph, a Feynman graph y, 
representing a specific contribution to Z, or rather to Z = Z/Zq. 

Z = £z y . (7) 

y 

The contribution Z r to the total Z from a arbitrary Feynman graph y defines the value of 
the graph, which is given by the following Feynman rules. 

1. Each vertex i is associated with a factor of v mp where m ; is its degree; 

2. Each edge is associated with a factor of b~ l ; 

3. The vertex and edge factors are multiplied together, and the result divided by the 
proper edge symmetry factor S y (and of course by AH). 

The edge symmetry factor S y is the usual one, accounting for teh order of the edge 
symmetry group, with a factor of n\ for each m-fold edge between distinct vertices, and 
a factor of 2 n nl for each vertex with n self-connections (where the 2" comes from the 
symmetry under flipping edges involved in self-connections). 

The above rules apply to labelled graphs, where vertices are considered distinguish- 
able, while edges are indistinguishable, as are their directions. In many cases one con- 
siders instead unlabelled graphs; the only difference is that then, also vertices are con- 
sidered indistinguishable. As a result, the 1/7Y! factor will be replaced by the correct 
vertex symmetry factor of the particlar graph considered, which together with the edge 
symmetry factor yields a total symmetry factor Sy. 

We will, however, focus on labelled graphs, that have distinguishable vertices. 



Saddlepoints and Degree Distributions 

Following refs. [13, 14], we now wish to interpret the Feynman graphs associated 
with a particular integral as random graphs, with a statistical weight proportional to the 
corresponding contribution to Z, as given by the Feynman rules above. 

Then it is interesting to investigate what kind of random graph ensembles this leads 
to. For a fixed N, it is obvious that the probability for a certain degree sequence, {mi}, 
is given by 

(8) 

Z N 

If we wish to focus on the degree of a single vertex, all we have to do is to sum over the 
N—l other degrees, to obtain the single vertex degree distribution 

v w /g~ fa ^ / v , y(*y v ~ 1 dx 

Pm ~ m\ fe- b * 2 / 2 V(x)"dx ( } 



having the generating function 



tj( \ _ y m Je-^^V(xz)V(xf^dx 

H ^ = l PmZ = Ie -^ V{x)Ndx < 10 > 

Now, assume for large N that Z/v can be evaluated in a saddlepoint approximation, 
obtained from the neighborhood of a real, positive saddle point x = xq, satisfying the 
saddle point condition corresponding to an extremal value of the "effective action", 
S N (x) = -bx 2 /2+N\ogV(x), 

S' N (x ) = bx -NV'(xo)/V(x ) =0 (11) 
To leading order, the saddle point approximation yields 

Z N oc e- x2 o/ 2 V» , (12) 

with Vq = V(xq), yielding 

Pm ^J§ * H lz) »^.. (13) 
m\Vo Vq 

Thus, the perturbation V(x) has an obvious relation to the generating function for the 
single-vertex degree distribution of Feynman graphs; this has been pointed out and 
explored to some extent, e.g. in ref. [13], where the perturbation was set to H from 
the start. Note, however, that any perturbation V with non-negative Taylor coefficients 
will do; by rescaling the integration variable, V can be recast into the form of H. 

We can also consider the joint distribution of the degrees of a small number of vertices. 
Thus, for the pair distribution, we obtain the nicely factorizable bivariate generating 
function, 

V ( x zi)V (x pzi) 

This factorization generalizes to larger sets, and indicates the approximate independence 
of the degree distributions of distinct vertices. 

To clear the view, let us change the integration variable to one rescaled by xq\ with 
y = x/xq and Vq = V(xq), the saddle point condition (11) yields b = NV'(xq)/xqVq, and 
we obtain, apart from an uninteresting constant, 

Z N = J exp (-^v 2 ) H(y) N dy (15) 

with an integrand designed to have an extremum at y = 1 ; m is defined as the expected 
average degree, as given by m = H'(\), while H(l) = Y*mPm = 1 of course is assumed. 
Note that the integrand has the form of a fixed factor taken to the Mh power, indicating 
that for large N, the saddle point approximation should be OK (provided the correct 
saddle point is used). 

In what follows, we will assume we have performed such a rescaling to a natural 
variable y, and that Z N is defined by eq. (15), with H(z) = YlmPmZ m taken to generate a 
desired degree distribution {p m }. 



H 2 (z h z 2 ) ~ " I2 ~H(zi)H(z 2 ) (14) 



Similarities and Differences to DRG 



Based on the expression eq. (15) for Z/v, let us compare the resulting ensemble of 
Random Feynman Graphs (RFG) and compare it to the corresponding DRG model with 
a degree distribution {p m } given by the one generated by H(x). 

Let us first note that Feynman graphs are based on Gaussian integrals, where the value 
(Af - 1 ) ! ! of the integral fx M e~ bx2 / 2 dx for an even M obviously stems from the number 
of distinct pairings of the M factors of x in x M . Thus, for a given degree sequence {mj 
yielding an even total stub count M = Y*i m u the distribution over compatible Feymnan 
graphs is obviously given by random stub pairing, just as in DRG. Thus, the restrictions 
to a fixed degree sequence of a DRG model and the associated RFG model are obviously 
identical, and so it is enough to compare the distributions over degree sequences, P {mi}. 
In DRG, this factorizes as Yli Pm P with the small modification that M must be even. 

For RFG, we obtain, with M = Y*i m u 



, 'fe-*T xZ H(x) N dx V "" v ' fe-f x2 H(x) N dx 

(16) 

which, apart from a factor depending only on M, is identical to the DRG expression. 
This implies that the restrictions to a fixed total stub count M of a DRG model and the 
associated RFG model are identical, and so it is enough to compare the distributions 
over total stub counts, Pm- 

The M-distribution is simplest expressed in terms of its generating function, U (z) = 
Y,mPmz M - For DRG, it is essentially given by U (z) ~ H(z) N . The trivial requirement of 
an even M yields the slightly modified expression 



H (z) N +H(-z ) N 



u ° R0 (z) = ii l< d7) 



corresponding to the distribution 



UPrm (18) 

{«,-} ' 

For RFG, on the other hand, the result is obviously given by P^ G = yielding the 
generating function 

fe-^ x2 H(x) N dx 
The resulting M-distribution can be written as 

oc (M - 1 ) ! ! (Nm) ~ M/2 P^ G , (20) 

and that is essentially the only difference between the models. 

Thus, upon clamping the total stub count M to an even number close to the expected 
value M = Nm, a pair of associated DRG and RFG models are identical [13, 14]. Note 



that the relative factor in eq. (20), (Nm)~ ' P^°, is stationary for the expected value 
M = Nm, which shows that the agreement for clamped M is not very sensitive to the 
precise clamped value, as long as it is not too far from the expected value M. 



Validity of the Saddle Point Approximation 

At the basis of the above discussion is the assumption that the saddle point con- 
tributions really dominates for large N, which is not obviously always the case. 
A necessary condition is that xq = 1 defines a local maximum of the effective action 
S eff (x) °c — mx 2 /2 + \ogH(x), yielding the rather weak requirement 

S'^(x= 1)« (m(m- 2))- (m) 2 < (m 2 ) - (m) 2 < 2 (m) (21) 

This does not appear to correspond to any obvious property of the associated graph 
ensemble. 

In case this condition is not fulfilled, the contributions from x ~ 1 will fail to dominate 
Z/v for large N. A possible scenario then is that a competing saddle point xq ^ I will 
take over; assume it is on the positive real axis. Then we can always define a rescaled 
integration variable x = x/xq, such that the proper saddle point is moved to x = 1. 
The resulting degree distribution {p m } will have a geometric factor as compared to the 
expected one, 

- _ H(xqx) _ x™p m 

H{x) ~ m^) * Pm ~ hW) • (22) 

Similarly, if we consider the full Z instead of merely the restriction to a fixed order 
N, we have, in the rescaled variables, 

Z - / e NoS ^dx, with (23) 



x 2 

s(x) = -m—+H(x) (24) 

where the parameter Nq is to be taken as a desired value of N. Then s'(\) =0, and 
we have a saddle point at x = 1, just as for Z^ with Af = A^o- This time, however, 
the requirement of a local maximum, s"(l) < 0, is stronger; furthermore, it has an 
obvious interpretation: it becomes the well-known condition of subcriticality, i.e. the 
non-existance of a giant component in the associated DRG model, 

-m + H"(x) < (m(m - 2)) <0 (25) 

It is also interesting, for a fixed large No, to compare the full Z to the restrictions Zm 
for different N in the neighborhood of Nq, and to the double restrictions for different 
M in the neighborhood of Mq = Nom. Of course, Z will almost always define a divergent 
integral, but we can assume it to be formally dominated by the saddle point value, and 
compare that to the contributions from different Af and M. We have 

Z = y e -^omx 2 /2+M)//W (k (26) 



Z N = ^Je- N ^l 2 H{x) N 6x (27) 

Z NM = %J^f MH ^ N fe-^Wdx (28) 

(29) 

where we assume Z to have a proper saddle point at x = 1, i.e. (m(m — 2)) < 0. This 
implies the saddle point value g^oO-m/ 2 f or z, as well as for Z/y with N = No. For Zm 
with a slightly different N, the value will to lowest order not change - the prefactor 
is obviously staionary, and due to the dominance of the saddle point x = 1, a slightly 
different power of H(l) = 1 makes no difference. 

Thus, Zn as a function of ./V will have an extremum at /V = No; furthermore, the 
corresponding saddle point value coincides will that of Z. In a similar way, Z NM as a 
function of M for fixed /V = A^o, will have an extremum for M = Nth, and the saddle 
point value coincides with that of Z#. This indicates that under certain circumstances, Z 
gets its most important contributions from N ~ Nq and M ~ NqM, as it should. 



RANDOM FEYNMAN GRAPHS - MULTIVARIATE CASE 

Now let us investigate instead the slightly more complicated case of a multivariate 
integral, representing a statistical mechanical partition function for a AT-dimensional 
variable x = {x a , a = 1 . . .K}. Thus, consider the ^-dimensional integral 

Z = [e s ^d K x = y>(x)+v(*) d * x = | e 4* T Bx e K(x) d K x (30) 

with an action S(x) in the form of an unperturbed quadratic action So(x) = 
— Y*ab x aBabXb/2 plus a perturbation V(x). We will restrict our attention to cases 
where (i) the inverse B _1 of the symmetric matrix B has no negative elements, and (ii) 
the perturbation V(x) has only non-negative multivariate Taylor coefficients, i.e. 

Y m a 

with v m > 

in terms of the multivariate power m = {m a ,a = 1 . . .K}, and m! = W a m a \. 

Feynman Graphs for a Multivariate Integral 

The partition function Z can be evaluated in a perturbative manner around the Gaus- 
sian approximation, in complete analogy to the univariate case, with the perturbative 
expansion organized in the form of a sum of contributions associated with Feynman 
graphs, as follows. 



V x 

'—' 111 



(31) 
(32) 



The first step is to expand e v ( x \ yielding the decomposition 

z = E§7' with ( 33 ) 
Z N = J e-? xTEx V(xfd K x (34) 

Then, each of the N factors V(x) is associated with a distinct vertex i, and Taylor- 
expanded to yield 

Z N = Z {m,} 5 wim (35) 

z w = fi v ^f^ xEx n- m ^ (36) 

For each of the factors x mi , its multivariate power m, is interpreted as the colored degree 
of vertex i; thus, each of its m ia factors x a is associated with a stub with color a, belonging 
to vertex i, which has the (plain) degree m,- = E a m ;a . 

It is obvious that the integral in eq. (36) depends on {m,} only via the total colored 
stub count M = £ ; m,. Again, as in the univariate case, the value can be seen as arising 
from all possible ways to pair the M = }2 a M a (which again must be even) distinct x- 
f actors. 

The pairings can be organized in equivalence classes, corresponding to distinct Feyn- 
man graphs, as usual characterized by a unique adjacency matrix {n*/}. Each element 
riij counts the number of stubs of vertex i that are paired with a stub from vertex j. The 

value of Z = Z/Zq, with Zq = J e~3 xTBx cbc, becomes the sum over contributions from 
the individual Feynman graphs y, 

/ Z'/y (37) 

y 

These counts can be further decomposed into the elements of a colored adjacency matrix 
{niajb}> w i m elements n ia j b counting the number of a-colored stubs of vertex i that are 
paired with a ^-colored stub of vertex j, such that the sum rules w (a = Y.jb n iajb and 
nij = Lab n iajb hold. 

Thus, the contribution to Z from each graph y - its value Z y - is due to a sum over 
pairings consistent with that graph, and is given by the following Feynman rules for the 
contribution. 

1. Each stub in y is assigned an independent color variable a. 

2. Each vertex i is associated with a factor v m/ , where m, = {m ja } is its colored degree; 

3. Each edge with stub colors (a, b) at its endpoints is associated with a factor of B^ 1 ; 

4. The vertex and edge factors are multiplied together, the result is summed over all 
stub colors, and the result divided by the proper edge symmetry factor S r (and of 
course by TV!). 



These rules are obvious generalizations of the univariate versions. 



Now we want to interpret these Feynman graphs as random graphs in a multivariate 
version of RFG, with a statistical weight for each distinct graph proportional to the 
corresponding graph value as given by the Feynman rules. 



Saddlepoints and Colored Degree Distributions 

We wish to investigate the properties of the resulting graph ensemble. To that end, let 
us investigate the degrees, or even better, the colored degrees, of the resulting graphs. 
For a fixed N, the distribution over colored degree sequences {m,} is obviously given 
by 

P{m ; -} = %^ (38) 

By summing over the colored degrees of all vertices except one, the colored degree 
distribution of a single vertex results, and is given by 

_ v m Je-^x m V(x) N - l d K x 

Pm— — : f-rr 

m! J e -2^V{x) N d K x 

As was the case for univariate integrals, we can hope to evaluate this with saddle point 
methods in the thermodynamic limit of N — > °°. For a point xo to be a saddle point of Z#, 
it must render the integrand of Z N extremal, corresponding to the saddle point condition 

VV 

= -Bx + ^V— =s> Bx V(x )=;VW(xo) (40) 
For the single-vertex colored degree distribution, the saddle point approximation yields 

Pm = (41) 
m!V(x ) 

with the generating function H(z) = L m p m z m given by 

where the ring "o" represents the componentwise multiplication of two vectors. 

Assuming the saddle point xo to have entirely positive components, xq m > 0, we can 
change integration variables to rescaled variables y, given by x = yoxo- Accordingly, we 
can define Vq = V(xo), in terms of which V(x) = V(yoxo) = VoH(y). We also choose 
to define a transformed matrix T from AT -1 = xo o B o xo. yielding T ab > 0. 

For Z N , this yields the expression 



je-^ TT ^H(y) N d K y (43) 



with an integrand constructed to have a saddle point aty = l = {l,...,l}, which implies 
the following constraint on the matrix T. 



T 1 = ^y Tm = l (44) 

where we have introduced the average colored degree, m = V//| y=1 , while normaliza- 
tion of the colored degree distribution with necessity implies H(\) = 1. 



Similarities and Differences to CDRG 

Now let us assume we are using the convenient coordinates x, with the saddle point 
asymptotically assumed to be at x = 1, in terms of which Z N is defined by eq. (43). 

From the above discussion, it is should be obvious that the multivariate RFG approach 
shows strong similarities to the corresponding CDRG model [10], which is defined in 
terms of a given graph size N, a color set [1, . . . ,K], a colored degree distribution {p m }, 
and a color preference matrix, T = {T ab }, with non-negative elements T ab , required to 
fulfill the constraint eq. (44), as follows. For each vertex, its colored degree is drawn 
independently from the given distribution. Then a weighted random pairing of the full 
set of M = Y. a M a colored stubs is performed, where each edge connecting colors a,b 
provides a weight factor oc T ab . 

For the random Feynman graphs, the distribution over colored degree sequences, for 
fixed N, becomes 

This can be compared to the corresponding CDRG result, which of course has the 
factorized structure 

^ = ( 46 ) 

i 

In analogy to the univariate case, we note that the associated RFG and CDRG ensembles 
agree on (1) the distribution over colored graphs conditional on a fixed colored degree 
sequence, and (2) the distribution over colored degree sequences conditional upon a 
fixed total colored stub count M. The former follows from the equivalence of the random 
pairing step involved in both ensembles, and the latter from a direct comparison of eqs. 
(45,46). 

Thus, the two ensembles only disagree on the distribution over the total colored stub 
count M; once it is clamped to a certain value, the two yield identical distributions over 
graphs, even before summing over colors. 

In agreement with the univariate case, the condition for the saddle point 
at x = 1 to define a local maximum of the integrand for the formal integral 
Z = fe N °(~ xTT lx /2+#(x)) d # x is e q U i va i ent to t h e condition of subcriticality of the 
associated CDRG model [10]. Also, Z is formally dominated by graphs with N = Nq 



and M = Nom, in the sense that the saddle point values of Z and the restrictions Zm and 
Z/vm all agree to leading order, while Z N as a function of ./V is stationary for N = Nq, as 
is Z Nq m as a function of M for M = Nq (m) . 

CONCLUSIONS 

A large class of ensembles of Feynman graphs for integrals representing statistical- 
mechanical partition functions have been investigated as models for random graphs, 
with the statistical weight of a graph taken to be proportional to its associated value as a 
contribution to the integral. 

For the case of univariate integrals, it was previously known [13] that the resulting 
random Feynman graph ensembles closely resemble those of random graphs with a 
given degree distribution, or DRG. In this article, it was found that for the case of 
multivariate integrals, the associated random Feynman graphs in a similar way resemble 
those of CDRG, a stub-colored extension of DRG. These findings are new, although such 
relations were conjectured in ref. [10]) 

For the restriction of the ensembles to a clamped value of the total (colored) stub 
count, the agreement was in fact shown to be exact, while the distribution over the total 
(colored) stub count was shown to differ, which leads to differences also in the resulting 
degree distributions. We speculate that this may be the effect of the random Feynman 
graph ensemble defining a kind of annealed approximation to the (C)DRG one, where 
the disorder as defined by the random (colored) degree sequence is considered quenched, 
i.e. randomly chosen from a given distribution, then clamped in the random pairing step. 
In the random Feynman graphs, as a contrast, both degrees and pairings fluctuate on an 
equal level. 

The establishment of close relations between different models opens up the stage 
for questions as to the relation between the critical phenomena in the graph models 
(such as the appearance of a giant component) and those in the associated statistical- 
mechanical models. Some relations of this type was found, but more work is needed to 
further illuminate this issue. 

Finally, we note that a straightforward generalization to random Feynman graph 
models for directed stub-colored graphs should be possible by considering suitably 
defined multidimensional integrals over a set of complex variables. 
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